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We report on a measurement of the Cabibbo-Kobayashi-Maskawa CP-violating phase 7 through a 
Dalitz analysis of neutral D decays to Ks1t~it + in the processes 

B=F _> d { * ) K+, D* -f Dtt°,D^. 

Using a sample of 227 million BB pairs collected by the BABAR detector, we measure the amplitude 
ratios r B = 0.12±0.08±0.03±0.04 and r* B = 0.17±0.10±0.03±0.03, the relative strong phases 8 B = 
(104±45 +H tit) and 8% = (-64 ± 41 +_\ 4 2 ± 15)° between the amplitudes A(B~ -> D^ K~) 
and A(B~ -> D W °A"~), and 7 = (70 ± 31 +_\l t\t)° ■ Tne first error is statistical, the second is 
the experimental systematic uncertainty and the third reflects the Dalitz model uncertainty. The 
results for the strong and weak phases have a two-fold ambiguity. 

PACS numbers: 13.25.Hw, ll.30.Er, 12.15.Hh, 13.25.Ft 



CP violation in the B meson system has been clearly 
established in recent years [1, 2]. Although these re- 
sults are in good agreement with Standard Model ex- 
pectations, other and more precise measurements of CP 
violation in B decays are needed to over-constrain the 
Cabibbo-Kobayashi-Maskawa (CKM) quark mixing ma- 
trix [3] and search for new physics effects. The angle 7 of 
the unitarity triangle [4] of the CKM matrix constitutes 
one of these crucial measurements. 

Various methods using B~ -> D(*)°K- [5] decays have 
been proposed to measure 7 [6-8]. Here, D° indicates ei- 
ther a D° or a D° meson and the symbol "(*)" refers to 
either a D or D* meson. All methods exploit the fact that 
a B~ can decay into a D^°R- (D^K') final state via 
b — > cms (b — ► ucs) transitions. These decay amplitudes 
interfere when the D° and D° decay into the same fi- 
nal state, which can lead to different B+ and B~ decay 
rates (direct CP violation). In this Letter we report on a 
measurement of 7 based on the analysis of the Dalitz dis- 
tribution of the three-body decay D° — > K^tt~tt+ [7, 8]. 
The primary advantage of this method is that it involves 
the entire resonant structure of the three-body decay, 
with interference of doubly Cabibbo-suppressed (DCS), 
Cabibbo-allowcd (CA), and CP eigenstate amplitudes, 
providing the sensitivity to 7. The analysis is based on an 
integrated luminosity of 205 fb" 1 recorded at the T(4S) 
resonance (corresponding to 227 million BB decays) and 
9.6 fb _1 collected at a center-of-mass (CM) energy 40 
MeV below with the BABAR detector [9] at the SLAC 
PEP-II e+e~ asymmetric-energy B Factory. 

The small CP asymmetry in D decays allowed by 
the present experimental limits [10] has a negligible ef- 
fect on this analysis. Thus, the B+ -> D^°K+, 
D*° -> L»°7r°,I) 7, D° -> K a s TT-TT+ decay chain ampli- 
tude A^\rn 2 _,m 2 _) can be written as 

A D (m 2 T ,m 2 ± ) + Kr^e^B+^ADim^ml) , 

where m 2 _ and mZy are the squared invariant masses 
of the K®n~ and K^tt+ combinations, respectively, 
and AD(m 2 -,rn 2 _) is the D° — > K®ir~Tr+ decay am- 
plitude. Here, r$ and 5^ arc the amplitude ra- 



tios and relative strong phases between the amplitudes 
A(B~ -> DW°K-) and A(B~ -> D^°K~). As a con- 
sequence of parity and angular momentum conserva- 
tion in the D*° decay, the factor k takes the value 
+1 for B~ -» D°K- and B~ -> D* (D°ir°)K- , and 
-1 for B~ -> D*°(D a "/)K- [11]. We first deter- 
mine ^4jj(m?_, m^_) through a Dalitz analysis of a high- 
statistics sample of tagged D° mesons from inclusive 
D*+ — > D Q n+ decays reconstructed in data. We then 
perform a simultaneous fit to the \A < "*\m 2 _, m 2 + )\ 2 and 
\A { *\m 2 _,m 2 + )\ 2 distributions for the B+ -> D^°K+ 

samples to determine the CP parameters r y B , S B , and 
7. We emphasize that in this analysis the Dalitz ampli- 
tude is only a means to extract the CP parameters. 

B~ candidates are formed by combining a mass- 
constrained £>W° candidate with a track identified as 
a kaon [9]. We accept K® — > tt + it~ candidates that have 
a two-pion invariant mass within 9 McV/c 2 of the K® 
mass [4] and a cosine of the angle between the line con- 
necting the D° and K® decay vertices and the mo- 
mentum (in the plane transverse to the beam) greater 
than 0.99. D° candidates are selected by requiring the 
K®ir~Tr+ invariant mass to be within 12 McV/c 2 of the 
D° mass [4]. The tt° candidates from D*° -> D°ir° are 
formed from pairs of photons with invariant mass in the 
range [115, 150] MeV/c 2 , and with photon energy greater 
than 30 MeV. Photon candidates from D*° -> L>°7 
are selected if their energy is greater than 100 MeV. 
D*° — > D tt (D°j) candidates are required to have a 
D*°-D° mass difference within 2.5 (10) MeV/c 2 of its 
nominal value [4] . 

The beam-energy substituted B mass ttzes [12] 
(Fig. 1) and the difference AE between the reconstructed 
energy of the B~ candidate and the beam energy in the 
e+er CM frame are used to identify signal B~ decays. 
We require m ES > 5.2 GeV/c 2 and \AE\ < 30 MeV. 
Since the background is dominated by random combina- 
tions of tracks arising from e+e~ — > qq, q = {u,d, s,c} 
(continuum) events, we require | cos 6^1 < 0.8, where 0? 
is the CM angle between the thrust axis of the B~ can- 
didate and that of the remaining particles in the event. 
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FIG. 1: Distributions of m ES for (a) S" -> D° 'K~ , (b) B~ -> 
D* {D tt°)K~, and (c) B~ -> D*°(D°-f)K~. The curves 
represent the fit projections for signal plus background (solid 
lines) and background (dotted lines). The peaking structure 



of the background is due to remaining B — > _D'*-"V events 



The reconstruction efficiencies (purities in the signal re- 
gion m ES > 5.272 GeV/c 2 ) are 18% (63%), 5.9% (86%), 
8.1% (52%) for the B~ -> D°K-, B~ -> D* a (D a n a )K^ , 
and £>~ — > D*°(D°j)K~ decay modes, respectively. The 
cross-feed among the different samples is negligible. 

The D° decay amplitude is determined from an un- 
binned maximum-likelihood Dalitz fit to a high-purity 
(97%) sample of 81496 D*+ -> D vr + decays recon- 
structed in 91.5 fb _1 of data (Fig. 2). We use the isobar 
formalism described in Ref. [13] to express ylc as a sum 
of two-body decay-matrix elements (subscript r) and a 
non-resonant (subscript NR) contribution, 

A D {m 2 _,m\) = E r a r e i4,r Ar{m 2 _,m 2 + ) + a N Re l0NR , 

where each term is parameterized with an amplitude 
a r and a phase <p r . The function A r (m 2 _ 1 m 2 Sr ) is the 
Lorentz-invariant expression for the matrix element of a 
D° meson decaying into K^tt~tt + through an intermedi- 
ate resonance r, parameterized as a function of the posi- 
tion in the Dalitz plane. 

Table I summarizes the values of a r and <p r obtained us- 
ing a model consisting of 16 two-body elements compris- 
ing 13 distinct resonances and accounting for efficiency 
variations across the Dalitz plane and the small back- 
ground contribution. For r = p(770), p(1450) we use the 
functional form suggested in Ref. [14], while the remain- 
ing resonances are parameterized by a spin-dependent 
relativistic Breit-Wigner distribution. For intermediate 
states with a K* , the regions of interference between DCS 
and CA decays are particularly sensitive to 7, and we 
include the DCS component when a significant contribu- 
tion is expected. In addition, we find that the inclusion of 
the scalar 7T7r resonances a and a' significantly improves 
the quality of the fit [15]. Since the two a resonances are 
not well established and are only introduced to improve 
the description of our data, the uncertainty on their exis- 
tence is considered in the systematic errors. We estimate 
the goodness of fit through a two-dimensional x 2 test and 
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FIG. 2: (a) The D° -» KIiy~-k + Dalitz distribution from 
D* + — > D°7r + events, and projections on (b) m 2 _, (c) m+, 
and (d) m^ +jr _. The curves are the fit projections. 



TABLE I: Amplitudes a r , phases cj> r and fit fractions ob- 
tained from the fit of the D° — > K^-k~-k + Dalitz distribu- 
tion from D* + — > D°tv + events. Errors are statistical only. 
Masses and widths of all resonances except a and a' are 
taken from [4]. The fit fraction is defined as the integral of 



:\A r (r 



, m 

of \A D (m 2 _ ,m\)\ 2 . The sum of fit fractions is 1.24. 



over the Dalitz plane divided by the integral 



Resonance 


Amplitude 


Phase (deg) 


Fit fraction 


A"*(892)" 


1.781 ±0.018 


131.0 ±0.8 


0.586 


KS (1430) " 


2.45 ±0.08 


-8.3 ±2.5 


0.083 


if 2 *(1430)- 


1.05 ±0.06 


-54.3 ±2.6 


0.027 


#-(1410)- 


0.52 ± 0.09 


154 ± 20 


0.004 


A"*(1680)" 


0.89 ±0.30 


- 139 ± 14 


0.003 


/T(892) + 


0.180 ±0.008 


-44.1 ±2.5 


0.006 


ir *(1430)+ 


0.37 ±0.07 


18 ±9 


0.002 


K* 2 (1430)+ 


0.075 ± 0.038 


- 104 ± 23 


0.000 


p(770) 


1 (fixed) 


(fixed) 


0.224 


w(782) 


0.0391 ± 0.0016 


115.3 ± 2.5 


0.006 


/o(980) 


0.482 ± 0.012 


-141.8 ±2.2 


0.061 


/o(1370) 


2.25 ±0.30 


113.2 ±3.7 


0.032 


/ 2 (1270) 


0.922 ± 0.041 


-21.3 ±3.1 


0.030 


p(1450) 


0.52 ±0.09 


38 ± 13 


0.002 


a 


1.36 ±0.05 


-177.9 ±2.7 


0.093 


a 


0.340 ± 0.026 


153.0 ±3.8 


0.013 


Non Resonant 


3.53 ± 0.44 


128 ±6 


0.073 



obtain \ 2 = 3824 for 3054 — 32 degrees of freedom. 

We simultaneously fit the B~ — > D^*'°K~ sam- 
ples using an unbinned extended maximum-likelihood 
fit to extract the CP-violating parameters along with 
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the signal and background yields. Three different back- 
ground components are considered: continuum events, 
B~ -> D^°n- and T(45) -> BB (other than B~ -> 
2}(*)°7r~) decays. In addition to rngs, the fit uses AE 
and a Fisher discriminant [12] to distinguish signal from 
B~ — > D^^n^ and continuum background, respectively. 
The log-likelihood is 



ln£ = -JI^ + X! 111 



E^(S)^ Dalitz (*) 



where £j = {toesj &E,!F}j and ry} = {m 2 _,m^_)j char- 
acterize the event j. Here, P c (0 an d 'P° alltz (?7) are the 
probability density functions (PDF's), and iV c the event 
yield for signal or background component c. For signal 
events, P^ ahtz (ff) is given by |.4^(f7)| 2 corrected by the 
efficiency variations. All PDF shape parameters used to 
describe signal, continuum and B~ — > Z)(*'°7r~ compo- 
nents are determined directly from B~ — > D^*'°K~ and 
7r signal, sideband regions, and off-peak 
data, and are fixed in the final fit for CP parameters and 
event yields. Only the ttiesj A23 and Dalitz PDF's for 
BB background events are determined from a detailed 
Monte Carlo simulation. B~ — » Z?(*)°7r~ candidates have 
been selected using criteria similar to those applied for 
B~ — > L)(*' K~ but requiring the bachelor pion not to 
be consistent with the kaon hypothesis. 

The CP fit yields 282 ± 20, 90 ± 11, and 44 ± 8 signal 
D°K-, D*°(D°ir°)K-, and D*°(D°j)K- candidates, 
respectively, consistent with expectations based on mea- 
sured branching fractions and efficiencies estimated from 
Monte Carlo simulation. The results for the CP-violating 
parameters = (x± , y± ), where and y± are 



defined as the real and imaginary parts of the complex 

±7 ) respectively, are summa- 



amplitude ratios r^±e*^ ~ 



(*) 

rized in Table II. Here, r B ± is the amplitude ratio be- 
tween the amplitudes b — ► u and b — > c, separately for 
B + and B~ . The only non-zero statistical correlations 
involving the CP parameters are for the pairs z_, z + , zl, 

and z^J_, which amount to 3%, 6%, —17%, and —27%, 

(*") 

respectively. The Zj_ variables are more suitable fit pa- 

rameters than r B , <$£ and 7 because they are better 
behaved near the origin, especially in low-statistics sam- 
ples. Figures 3(a,b) show the one- and two-standard de- 
viation confidence-level contours (statistical only) in the 
zW planes for D°K~ and D*°K~ , and separately for B~ 
and B + . The separation between the B~ and B + regions 
in these planes is an indication of direct CP violation. 

The largest single contribution to the systematic un- 
certainties in the CP parameters comes from the choice of 
the Dalitz model used to describe the D° — > K®tt~it + de- 
cay amplitudes. To evaluate this uncertainty we use the 
nominal Dalitz model (Table I) to generate large samples 
of pseudo-experiments. We then compare experiment by 
experiment the values of Zj_ obtained from fits using the 



TABLE II: CP-violating parameters z^* obtained from the 
CP fit to the B~ -> D^ )0 K~ samples. The first error is sta- 
tistical, the second is the experimental systematic uncertainty 
and the third reflects the Dalitz model uncertainty. 





x ± 




z_ 


0.08 ±0.07 ±0.03 ±0.02 


0.06 ±0.09 ±0.04 ±0.04 


Z-\- 


-0.13 ±0.07 ±0.03 ±0.03 


0.02 ±0.08 ±0.02 ±0.02 


z* 


-0.13 ±0.09 ±0.03 ±0.02 


-0.14 ±0.11 ±0.02 ±0.03 


4 


0.14 ±0.09 ±0.03 ±0.03 


0.01 ±0.12 ±0.04 ±0.06 
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FIG. 3: Contours at 39.3% (dark) and 86.5% (light) confi- 
dence level (statistical only) in the z'*' planes for (a) D°K~ 
and (b) D*°K~ , separately for B~ (thick and solid lines) and 
B + (thin and dotted lines). Projections in the r^' — 7 planes 
of the five-dimensional one- (dark) and two- (light) standard 
deviation regions, for (c) D°K~ and (d) D*°K~. 



nominal model and a set of alternative models. We find 
that removing different combinations of K* and p reso- 
nances (with low fit fractions) , or changing the functional 
form of the resonance shapes, has little effect on the total 
X 2 of the fit, or on the values of z±\ However, models 
where one or both of the a resonances are removed lead 

to a significant increase in the x 2 of the fit. We use the 

(*) 

average variations of z± corresponding to this second set 
of alternative models as the systematic uncertainty due 
to imperfect knowledge of Ad- 

The experimental systematic uncertainties include the 
errors on the toes> AE, and T PDF parameters for sig- 
nal and background, the uncertainties in the knowledge 
of the Dalitz distribution of background events, the effi- 
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ciency variations across the Dalitz plane, and the uncer- 
tainty in the fraction of events with a real D° produced 
in a back-to-back configuration with a negatively-charged 
kaon. Less significant systematic uncertainties originate 
from the imprecise knowledge of the fraction of real D°'s, 
the invariant mass resolution, and the statistical errors 
in the Dalitz amplitudes and phases from the fit to the 
tagged D° sample. The possible effect of CP violation in 
B~ — » £)(*)°7r~ decays and BB background was found 
to be negligible. 

A frequentist (Neyman) construction of the confidence 
regions of p = (rg, r* B , 5b, 5 b , 7) based on the constraints 
on has been adopted [4]. Using a large number of 
pseudo-experiments corresponding to the nominal CP fit 
model but with many different values of the CP fit pa- 
rameters, we construct an analytical (Gaussian) param- 
eterization of the PDF of z± as a function of p. For a 
given p, the five-dimensional confidence level C = 1 — a 
is calculated by integrating over all points in the fit pa- 
rameter space closer (larger PDF) to p than the fitted 
data values. The one- (two-) standard deviation region 
of the CP parameters is defined as the set of p values for 
which a is smaller than 3.7% (45.1%). 

Figures 3(c,d) show the two-dimensional projections 

(*) 

in the r B — 7 planes, including systematic uncertain- 
ties, for D°K- and D*°K~. The figures show that 
this Dalitz analysis has a two- fold ambiguity, (7, 5 B ) — > 

(7 + 180°, 5 { B ] + 180°). The significance of direct CP 
violation, obtained by evaluating C for the most proba- 
ble CP conserving point, corresponds to 1.6, 2.1, and 2.4 
standard deviations, for D°K~ and D* a K~, and their 
combination, respectively. Similar results are obtained 

using a Bayesian technique with uniform a priori proba- 
ta) (*) 

bility distributions for r B , 5 B and 7. 

In summary, we have measured the direct CP- violating 
parameters in B~ — > D^°K~ using a Dalitz analysis of 
D° -> K°tt-tt+ decays, obtaining r B = 0.12 ± 0.08 ± 
0.03±0.04 [0,0.28], r* B = 0.17±0.10±0.03±0.03 [0,0.35], 
S B = (104 ±45 tl\ S% = (-64 ±41 ^±15)°, 

and 7 = (70 ±31 t\l ±\$)° [12°, 137°]. The first error 
is statistical, the second is the experimental systematic 
uncertainty and the third reflects the Dalitz model un- 
certainty. The values inside square brackets indicate the 
two-standard deviation intervals. The results for 7 from 
B~ D°K- and B~ -> D*°R- alone arc (70 ± 38)° 
and (71 ±35)°, respectively (statistical errors only). The 
constraint on 7 is consistent with that reported by the 

Belle Collaboration [8] , which has a slightly better statis- 

(*) 

tical precision since our r B constraint favors smaller val- 
ues. We are grateful for the excellent luminosity and 
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